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Abstract 

Background, aim, and scope In this communication we 
reexamine the use of various least squares techniques 
(namely ordinary least squares (OLS), data least squares 
(DLS), and total least squares (TLS)) for the purposes of 
allocation as proposed in an earlier article in this journal 
(Marvuglia et al. (Int. J. Life Cycle Assess. 15:1020-1040, 
2010)). These methods are placed within the context of 
traditional methods of partitioning allocation. An equiva¬ 
lence between least squares techniques and traditional 
partitioning is noted and demonstrated on previously pub¬ 
lished brick production data. 

Materials and methods A short summary of the relevant least 
squares techniques is provided followed by a description of 
the problem of inventory calculation for the case of more 
products than processes. This is presented in terms of non¬ 
unique solutions to underdetermined systems of linear equa¬ 
tions with intensity matrix as unknown. 

Results and discussion We provide another analysis of the 
Sicilian brick production case study. Upon reexamination of 
the brick data, a number of disparities in the published inven¬ 
tories for brick can be more fully explained in terms of (1) data 
quality, (2) additional assumptions made to extract a solution 
of an underdetermined system of equations, and (3) discrep¬ 
ancy vectors. 
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Conclusions and outlook Like other types of partitioning, 
inventories produced by least squares techniques correspond 
to one of the non-unique solutions of an underdetermined 
system of linear equations. Based on this insight, we advise 
against the use of least squares techniques as a black-box 
approach to the allocation problem and conclude that the 
recommendation of TLS on the basis of its asymptotic prop¬ 
erties is not theoretically justified. 
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Partitioning • Regression 

1 Introduction 

A source of subjectivity in the inventory phase, allocation has 
been called one of the Achilles heels in the practice of life cycle 
assessment (Marvuglia et al. 2010, p. 1036). The life cycle 
inventory literature contains a variety of proposals for carrying 
out allocation, and Chembini and Stromman (2011) noted that 
expert preferences appear to be evenly distributed across the 
range of methods. Wardenaar et al. (2012) observed that a 
number of existing environmental protocols actually recom¬ 
mend different methods of allocation during the inventory 
phase and debated the implications for policy applications. 

Recently, Marvuglia et al. (2010, 2012) have proposed the 
use of several least squares techniques for the purposes of 
allocation and demonstrated their calculation on data 
representing a brick production network in Sicily. One part 
of the appeal of the least squares techniques seemed to be that 
the scale factors were unique or “optimal” approximate solu¬ 
tions of an overdetermined system of linear equations with 
respect to some loss or objective function, perhaps eliminating 
some subjective elements in the inventory calculation. They 
ultimately recommended a total least squares (TLS) approach 
on the basis of asymptotic properties (i.e., strong consistency) 
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of the TLS estimator for scale factors. Regarding the relation¬ 
ship between least squares methods and typical partitioning, 
they opined: 

“One could argue that solving the inventory problem 
directly in its rectangular form (without making the 
system matrix square) intrinsically leaves unallocated 
the single functions of a multi-functional process with¬ 
out actually solving the problem of allocation, but just 
hiding it. However, in our view, making resort to allo¬ 
cation has no other usefulness than allowing the trans¬ 
formation of a rectangular system into a square one”. 
(Marvuglia et al. 2010, p. 1024) 

In the next section, we summarize earlier work on the 
subject of partitioning of a matrix-based life cycle inventory 
(Cruze et al. (2013b)) which presented partitioning allocation 
as a particular solution of an underdetermined system of linear 
equations. Allocation is an attempt to state life cycle intensi¬ 
ties of one of several products given fewer observations than 
unknowns. From this viewpoint, inventories derived by least 
squares techniques are not unique solutions. Rather, they are a 
just a few among infinitely many solutions to an 
underdetermined system of linear equations. Furthermore, 
least squares techniques are not something done in lieu of 
partitioning. They are equivalent to partitioning, sometimes 
with weights that don’t make practical sense. These claims are 
illustrated through reexamination of the original brick data in 
Marvuglia et al. (2010). We conclude in Section 4 that least 
squares techniques are not a black box solution to the alloca¬ 
tion problem. 


2 Inventory calculation and allocation 

In our notation throughout this paper, a bold capital letter will 
denote a matrix. Bold lowercase letters denote column vec¬ 
tors, and an italic variable denotes a scalar value. Let n denote 
the number of processes,/; denote the number ofproducts, and 
q denote the number of distinct interventions included in the 
study. Partitioning in matrix LCI and the least squares tech¬ 
niques deals specifically with the case n<p , therefore at least 
one multi-output process must be present. 

A fundamental relation of life cycle inventory describes the 
net outputs of a system f pXl in terms of the scale i at which 
each process (column of the technology matrix A pXn ) must 
operate: 

As=f. (1) 

One may wish to know at what scale each process must be 
operated to provide some quantity of the j th product, the 
functional unit fj. Provided that the technology matrix is 


square with p=n and A invertible, the corresponding list of 
scale factors Sj may be found straightforwardly by inversion: 

Sj=A~ i f J . (2) 

If total interventions at each process are in direct proportion 
to the scale of operation, the life cycle inventory of the / h 
product, gj is given by the inner product of the interventions 
matrix B ?x „ with the vector of scale factors as follows: 

Sj = Bsj (3) 


For the invertible technology matrix, Heijungs and Suh 
(2002) define the intensity matrix as 

A = BAT X . (4) 

The inventory for one unit of the / h product, gj , is the f h 
column of A qxp . Trivially, the intensity matrix A must satisfy 

A T A r = B t . (5) 

When there are more products or functions present than 
processes, Eq. 1 is an overdetermined system of equations, 
and it may not have an exact solution. This marks the point of 
departure for much of the allocation literature. Many methods 
of partitioning have been developed in which multifunctional 
processes are broken into virtual, single-function processes to 
obtain an invertible technology matrix A, (of dimension p x /;), 
and shares of inputs and interventions are assigned in a 
partitioned interventions matrix B t (q x p) in proportion to 
some characteristic of the joint products (e.g., mass, economic 
value, etc.) so that the final inventory g * may be calculated as 

g* = ( 6 ) 


2.1 Three least squares methods 

Marvuglia et al. sought a solution of the overdetermined Eq. 1 
without having to explicitly define some new technology 
matrix A t and interventions matrix B t . They proposed ordi¬ 
nary least squares (OLS) allocation, data least squares (DLS) 
allocation and total least squares (TLS) allocation as three 
approaches to obtain an approximate solution (vector of scale 
factors) to Eq. 1 . All three types of least squares problem can 
be written compactly as variants of the TLS minimization 
problem (DeGroat and Dowling (1993); Van Huffel and 
Vandewalle (1991); Van Huffel (1997)): 
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min S:A A,Af\\[AA\Af]\\ F . 

subject to (A + AA)s = f + Af. ' 

The OLS result may be obtained by setting A/t=0, while 
the DLS result is obtained by setting A/=0 in Eq. 7 and then 
solving the programming problem (DeGroat and Dowling 
(1993)). In words, these programs seek a matrix of the 
smallest corrections to the technology matrix and functional 
unit, denoted [AA\Af\, in the sense that its Frobenius norm, 
denoted || , || Fj is minimal. Under these minimal changes, the 
input-output relationships of the perturbed system shown in 
the constraints of Eq. 7 may be satisfied exactly. Given scale 
factors as unique minimizers of Eq. 7, the inventory vector 
follows the computation shown in Eq. 3 since the scale factors 
are obtained in the dimension of processes. Since the three 
minimization problems differ from one another in terms of 
their objective functions and constraints, the scale factors and 
resulting inventories produced by the three least squares 
methods need not be identical. 

2.2 Allocation and the underdetermined system of linear 
equations 

When the number of products in the system is greater than the 
number of processes, at least one process is multifunctional, 
and Eq. 5 represents an underdetermined system of equations, 
with the number of equations equal to the number of processes 
contained in the technology matrix. Without loss of generality, 
one can consider just one type of intervention at a time so that 

A T n = bi, (8) 


where the vectors rj and bf are the / th rows of the intensity 
matrix A and interventions matrix B , respectively. 

When finding a solution via partitioning, the system of 
equations 

A T t r = b * (9) 

is solved by inversion to obtain a unique intensity vector, say 
r„, and the linear combination f r t is interpreted as a compo¬ 
nent of the inventory related to the functional unit. While the 
solution to Eq. 9 is unique, Eq. 8 has infinitely many exact 
solutions. The solutions to Eq. 9 and Eq. 8 are closely related, 
the unique solution of the former being just one particular 
solution of the latter. 

More in-depth study of the equivalence of inventories 
derived from those systems of equations is given in Craze 
et al. (2013b), but we wish to summarize the following here: 


1. A partitioning allocation as in Eq. 9 obtains one solution 
of Eq. 8 made under additional equality constraints called 
side conditions, 

2. Obtaining the solution under side conditions is equivalent 
to solving Eq. 8 by some generalized inverse of the 
technology matrix, A T The meaning of these two points 
may be made more concrete by study of OLS allocation 
since the generalized inverse in point 2 is explicitly 
known. The Moore-Penrose pseudoinverse, (A')' is used 
to solve Eq. 8 so that 

r OL S = (A r ) + b. (10) 

Taken together with point 1, there is an equivalent 
partitioning of the technology and intervention matrices so 
that r Q LS can be obtained from the operations in Eq. 9. 

Similarly, equivalent partitionings can be found which 
contain DLS and TLS intensities for the target product in the 
resulting intensity vectors. The upshot is that the least squares 
techniques pursued by Marvuglia et al. are not something 
done in lieu of partitioning; they are each equivalent to 
partitioning. Moreover, they are not unique solutions in the 
sense of Eq. 8 despite the fact that the vectors of scale factors 
may be optimal in the sense of the programming problem 
stated in Eq. 7. There is no reason to expect that a least squares 
solution of Eq. 7 produces an inventory that is more physically 
meaningful, interpretable or actionable than more standard 
types of partitioning. Indeed, the debate as to which method 
of allocation is “right” persists because in practice, different 
allocations entail specific constraints on some parameters of 
the intensity vector in Eq. 8, and yet the outcome is interpreted 
without regard to those constraints. 


3 The brick production study revisited 

The raw, unallocated data for the Sicilian brick production 
system may be found in Marvuglia et al. (2010). The original 
process technology matrix consisted of 8 processes which 
exchanged 14 economic flows. Six multifunctional processes 
were present, with recycled inerts being an additional joint 
output of the brick production process. In addition, 21 types of 
environmental interventions were considered; these included 
several types of natural resources consumed as well as emis¬ 
sions to air and water, and three wastes. 

Marvuglia et al. compared five types of allocations. Three 
inventories were based on the selection of an approximate 
solution of Eq. 1 by OLS, DLS, and TLS regression-based 
techniques applied to the 14><8 system. The other two 
partitioning-based allocations were derived from technology 
matrices with dimensions of 13 products by 13 virtual pro¬ 
cesses. One of these allocations was based on physical con¬ 
siderations while the other was a partitioning allocation based 
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on the economic valuation of the joint products. The reduction 
in number of products in these traditional methods of alloca¬ 
tion came from a partial substitution of recycled inerts for 
sand. 

The disparities across allocations were significant. A few of 
the many disparities shown in the original tables (Marvuglia 
et al. 2010, Table 11) are summarized as follows: 

• The life cycle releases of emissions and wastes, and 
uses of resources assigned to 1 t of bricks by the OLS 
allocation are all minuscule in comparison to any of 
the other allocation procedures, the largest releases 
and consumption being measured on the order of 
joules or grams as opposed to tens or hundreds of 
megajoules or kilograms. We note that this is due to 
the choice of the Moore-Penrose inverse to solve 
Eq. 8, which produces the shortest length solution 
among all possible solutions of an underdetemined 
system of linear equations (Golub and Kahan 1965, 
p. 205). 

• By convention, coal consumption (a resource input 
from outside the system boundary) should be negative, 
yet all five methods produce positive results spanning 
ten orders of magnitude, with OLS allocation on the 
order of 5 J, physical allocation at 856 megajoules and 
economic allocation on the order of 12,200 
megajoules ! 

• Water requirements are of the same order of magni¬ 
tude under physical and economic allocations, but of 
opposite signs, as are life cycle emissions of CO 2 , CO, 
and non-methane volatile organic compounds 
(. NMVOC ). 

• Marvuglia, Cellura and Heijungs observed similarities 
between physical allocation and DLS and TLS allocations 
for a select few emissions, wastes, and resources and 
concluded that this lent some credence to the TLS tech¬ 
nique, but even the two traditional allocations agreed 
closely on emissions of phosphorus, nitrogen, and adsorb- 
able organohalides ( AOX ). 

It is important to understand these disparities with regard to 
possible issues of data quality, the additional conditions im¬ 
posed by partitioning and least squares allocation, and the 
provided discrepancy vectors. 

3.1 Inventory disparities and data quality 

The brick production process required inputs of virtually all 
other products in the system. The network featured seven 
upstream processes: production of crude oil, oil-derivatives, 
natural gas, electricity and heat, clay, sand and gravel, and 
biomass. These processes were the unique producers of all the 
necessary inputs within the system boundary; however, they 


do not satisfy system-wide needs operating together at posi¬ 
tive scale in the unallocated matrix. This relates specifically to 
energy inputs originating in the heat and electricity, natural 
gas, and oil derivative processes. If the assumption of linearity 
is to be believed, then these processes cannot be operated 
together at positive scale to satisfy their own needs, much less 
to satisfy needs of other processes in the network. The pro¬ 
cesses upstream of brick are the unique producers of their 
respective materials or energy inputs within the system bound¬ 
ary, but the unallocated electricity and heat, oil derivatives and 
natural gas processes cannot be rescaled to operate together in 
the forward direction as follows: 

1. The electricity/heat process requires 4.27 MJ natural gas 
to produce its joint outputs; the natural gas process in 
Table 1 has been rescaled accordingly. 

2. The natural gas process requires 154 MJ of diesel oil to 
operate at that level, the oil derivatives process shown in 
Table 1 has rescaled accordingly. 

3. The combined needs of the oil derivatives and natural gas 
processes exceeds the heat and electricity resources avail¬ 
able, as indicated by the negative values in the list of the 
net outputs for these three processes. 

The negative energy values indicate a shortfall of electricity 
and heat inputs going to the oil derivatives and natural gas 
processes. An even larger shortfall occurs when one attempts 
to couple the five remaining processes with these three. This 
suggests that: 

• the provided unit processes may be poorly representative 
of actual processes involved in the network, or 

• that some important energy processes may have been 
overlooked within the system boundary, or 

• if the measurements on each process are accurate, the 
phenomenon is poorly described by a linear model. 

This apparent shortfall of energy inputs could arise 
simply if data for one or more of the obtained processes 
were poor representatives of actual processes used in the 
Sicilian network or if some other gross error was present. 
However, it has far reaching implications for all of the 
calculated inventories. 

We note that several of the originally tabulated inventories 
(Marvuglia et al. 2010, Table 11) failed to maintain sign 
convention. As discussed above, all five inventories showed 
that the natural resource coal was a positive quantity, indicat¬ 
ing that it had been produced or otherwise its use was some¬ 
how avoided on net by the system. The possibility of sign 
change in the final inventory can be attributed to the existence 
of one or more negative scale factors. The originally published 
scale factors for each method of allocation attempted in 
(Marvuglia et al. 2010, Table 10) show a mix of positive and 
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Table 1 A physical inconsisten¬ 
cy in a submatrix of the original 
technology matrix. The result is a 
system wide shortfall of energy 
inputs throughout the network 


Processes 



Net output 


Electricity/heat 

Oil derivatives 

Natural gas 


Products Electricity (MJ) 

1 

-4.93E+00 

0 

-3.93E+00 


Heat (MJ) 

2.48E+00 

-1.74E+02 

-2.08E-01 

-1.72E+02 


Diesel oil (MJ) 

0 

1.54E+02 

-1.54E+02 

0 


Fuel oil (MJ) 

0 

1.54E+02 

0 

1.54E+02 


Natural gas (MJ) 

-4.27E+00 

0 

4.27E+00 

0 


negative signs. We note that in allfive cases, the negative signs 
correspond to processes (or virtual processes for physical and 
economic allocation) which relate to heat, electricity, natural 
gas, and the oil derivatives. 

In system boundary expansion, negative scale factors 
are sometimes interpreted as an avoided activity or a 
volume of activity that has been displaced by other pro¬ 
duction. In this particular case, the processes associated 
with negative scale factors are not additional processes 
that might be offset by production, but rather the unique 
producers of their respective energy types within the sys¬ 
tem boundary. The existence of negative scale factors in 
this example runs counter to a claim made in Heijungs 
and Tan (2010) that an inverse-positive technology matrix 
(one in which A,> 0 holds elementwise) could be obtain¬ 
ed by substitution and partitioning assumptions, and that 
it would preclude the use of structural path analysis 
techniques a la Stromman et al. (2009), results on fuzzy 
methods for uncertainty propagation (Tan (2008); 
Heijungs and Tan (2010); Cruze et al. (2013a)), and the 
Monte-Carlo algorithm of Peters (2007). 

These distinctions are illustrated in Figs. 1 and 2; 
Tables 4 and 5 in Section 1.2 of the supporting informa¬ 
tion were used to derive these figures. Each of these 
figures represents activity within the rescaled virtual 
network for physical allocation (Fig. 1) and a slight 
variation of the physical allocation which sets Wdiesei 
O ii=0.01 while keeping all other four allocation weights 
fixed (Fig. 2). Selected interventions at each virtual 
process are also shown; the reported life cycle intensities 
are the sums of these interventions at all 13 virtual 
process. Line thickness indicates the order of magnitude 
of each flow. The virtual processes shown in red in 
Fig. 1 are associated with a negative scale factor. A 
literal interpretation of the physical allocation would 
suggest that not only is coal a net product of the system, 
some virtual processes make multiple products, and the¬ 
se products originate in unexpected places, e.g., the 
diesel virtual process produces heat, electricity and 
crude oil, and seems to draw diesel oil from the natural 
gas process. Directed graphs for each of the five 


different inventories presented in Marvuglia et al. 
(2010) would be similar, and each precludes the use of 
some of the aforementioned analytical techniques for 
error propagation and inventory hybridization. By con¬ 
trast, an inverse-positive outcome like that shown in 
Fig. 2 maintains sign of all inventories, and each virtual 
process serves just a single function as expected. 

Using MATLAB’s powerful symbolic features (see 
Section 1.2 of the supporting information for more de¬ 
tails), we determined that inverse-positivity, a property 
which would at least preserve the signs of all final inven¬ 
tories, depends exclusively on the pair of weights 
assigned to electricity and diesel oil, denoted (w e iectricity, 
^diesel oii)> in a conventional partitioning of the brick data. 
Keeping weights assigned to three other virtual processes 
(n"'white clay? n ; s aiK |, w 0 ij vc ca kc) fixed at theii physical allo¬ 
cation values, the projection shown in Fig. 3 illustrates 
that as the threshold is approached from either side, the 
life cycle intensity of coal can be driven to ±oo since the 
allocated technology matrix approaches singularity for 
these allocations. Only (w e iectricity, n’ diese i oil) pairs on the 
side of the threshold corresponding to small values of 
^diesel oil produce inventories of coal which conform to 
sign assumptions. 

As we discuss in the next section, inventories calculat¬ 
ed by several methods of allocation, even inverse-positive 
outcomes, may differ anyway as each solution of Eq. 8 is 
derived under a different set of additional constraints. In 
this example, however, data quality may have added to 
the disparities among inventories since the inventories are 
not even bounded and may change sign. This explains 
why the economic and physical allocations often differed 
by orders of magnitude in the originally tabulated inven¬ 
tories. One key lesson this example provides is that none 
of the presented partitionings, least squares allocations 
included, can overcome a data quality issue. If it is known 
that all elements of the published technology matrix are 
accurate so that it is indeed a legitimate technology ma¬ 
trix, then this should lead us to question assumptions of 
linearity, in which case existing partitioning and least 
squares methods may not even be pertinent. 
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-447.82 kg 


15.98 kg 

^ | -32.02 MJ 


V 



ture of the network is changed so that after partitioning some processes 


3.2 Additional assumptions and least squares equivalent 
partitionings of the brick data 

Although derived from believable information, it is un¬ 
derstood that as a representation of any actual production 
structure, the technology matrix after allocation is an 
imaginary creation. Joint products of the same process 
may not be produced independently, and inputs, re¬ 
sources, emissions and wastes assigned to the virtual, 
single-function processes vary seemingly arbitrarily from 


one choice of partitioning to the next. Marvuglia et al. had 
hoped that some of the more subjective elements of allo¬ 
cation might be eliminated through the use of various 
least squares optimality criteria. 

Two comments are worth making here: 

• First, since allocation is simply a matter of finding one 
particular solution to Eq. 8, inventories reached by tradi¬ 
tional procedures or by least squares techniques are both 
made under a set of equality constraints on non- 
identifiable intensities. Since the data contain no 
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102.93 kg 


-0.01 MJ 
41.11 kg-^ 1 


-0.88 MJ 
49.00 kg^i 


MJ Coal 


^■■1. 1538 kg 

^ kg C02 


•0.01 MJ 
► 45.60 kg 



--3.00 MJ 

► 2.93 kg 


I .►'18.93 kg 


-0.17 MJ 

603.43 kg 


/ 

0.00 MJ 


Fig. 2 Rescaled virtual network for an inverse-positive allocation. In this case, each virtual process serves a single function and is the unique producer of 
its respective energy source or material 


information about these constraints, they can be 
established arbitrarily with far-reaching implications for 
the interpretation of the inventory vector. 

• As motivated in Section 2.2 the least squares procedures 
do not necessitate the actual construction of square ma¬ 
trices, but they are equivalent to doing so, sometimes with 
weights that correspond to indefensible allocations, allo¬ 
cations with corresponding weights that sum to 1, but fall 
outside the interval [0, 1], 


Both of these points are further explicated below. 

3.2.1 Side conditions 

For the brick data, the corresponding Eq. 8 is a system of 8 
equations in the 14 unknown life cycle intensities for each 
product, none of which are uniquely identified in the data at 
hand. The derivation of these additional conditions, or side 
conditions, from the allocated technology matrix is discussed 
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Fig. 3 Life cycle intensity of coal associated with 1 t of brick as a and w 0 | ive ca k e = 0.8 in keeping with allocation weights for physical 

function of allocation weights Webtricity and Wdiesei oil- In this figure, the allocation 

remaining allocation weights have been fixed at w white o iay = w san d = 0.5 


in more detail in Cruze et al. (2013b), but for illustrative 
purposes, we present one set of additional constraints imposed 
by physical allocation to obtain an intensity vector with re¬ 
spect to coal consumption: 

t inerts — -85r san( ] (11) 

^electricity - 2-48rheat = 2.562r na tu ra l gas - 0.0009 (12) 

t’white clay ^red clay (12) 

t sand — ^gravel (14) 

t’olivc cake 1" straw — 0.2436k{]icsc] 0 j| 0. 0008 (15) 

^fueloil — ^diesel oil- (1 6) 


These side conditions represent several assumed rela¬ 
tionships about life cycle coal consumption among the 
various products included within the network boundaries. 
To reduce the dimension of products, Marvuglia et al. 
made an assumption that only 85 % of recycle inerts could 
be substituted for sand. This is also a statement about the 
relative impacts of the two materials (Eq. 11). Physical 
allocation suggests that coal intensities between some pairs 
of joint products are indistinguishable, namely that the 
effects per kilogram of either clay are the same (Eq. 13), 
inputs of coal into sand and gravel production are identical 
(Eq. 14), and the life cycle consumption of coal by both oil 
derivatives is the same (Eq. 16). In Eqs. 12 and 15, the 
relative magnitudes of life cycle coal use also depend on 
the inputs into those two processes. These equations are 
merely one possible collection of six additional assump¬ 
tions under which all 14 product intensities may be obtain¬ 
ed. Table 6 provided in the supporting information de¬ 
scribes the entire family of allocation equivalent con¬ 
straints in terms of 85 % recovery rate of inerts, the five 
allocation weights for other multifunctional processes, and 
any type of 21 given environmental interventions. 

For a comparison, Table 2 shows intensity vectors 
r phys ., /"econ., and r OLS • (These are the first rows of their 
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Table 2 Individual product inventories of life cycle coal (MJ) for brick 
production network dictated under the side conditions of physical, eco¬ 
nomic, and OLS allocations. Negative sign indicates net use by a product, 
positive sign suggests net avoidance 


Life cycle coal consumption (MJ) 

Physical 

Economic 

OLS 

^electricity 

4.08E-01 

3.22E+00 

-8.80E-03 

^heat 

4.11E-02 

1.21E+00 

1.60E-02 

^white clay 

-1.10E-03 

9.52E-03 

-3.85E-03 

^red clay 

-1.10E-03 

8.65E-03 

-1.38E-03 

f recycled inerts 

-2.61E-03 

2.24E-02 

3.28E-04 

^"sand 

-3.07E-03 

2.64E-02 

-7.24E-03 

^gravel 

-3.07E-03 

1.93E-02 

-6.51E-03 

^"olive cake 

-6.08E-05 

1.02E-02 

-9.80E-04 

^straw 

-1.52E-05 

4.12E-03 

-3.44E-04 

^"crude oil 

2.02E-03 

5.91E-02 

7.79E-04 

^diesel oil 

3.26E-03 

3.87E-02 

1.88E-04 

^*fuel oil 

3.26E-03 

1.52E+00 

-3.83E-02 

^natural gas 

1.20E-01 

1.46E+00 

7.58E-03 

^bricks 

8.56E+02 

1.22E+04 

4.75E-06 


respective intensity matrices.) Since coal is listed as a 
resource used by most of the given processes when oper¬ 
ated at positive scale, the final intensities for coal are 
expected to be negative indicating resource consumption. 
Only a few of the products match this expectation under 
physical and OLS allocation, and none of the products 
match expectation under economic allocation. The posi¬ 
tive intensities seem to suggest that these products some¬ 
how create or contribute to the avoided use of coal else¬ 
where. Furthermore, under the economic allocation, all of 
the products system wide appear to produce (or avoid the 
use of) net amounts of coal despite the fact that it is a 
direct input into seven of the eight processes. Finally, the 
combined life cycle consumption of coal by straw and 
olive cake (r st raw +r oiive cake) under OLS allocation is 
smaller than the direct consumption of coal into the whole 
biomass process (Marvuglia et al. 2010, Table 2). 
Furthermore, the biomass depends on the use of diesel 
oil as an input in the agricultural phase. Both processes 
must be operated to provide necessary inputs of biomass 
for the brick process, and yet, under OLS assumptions, 
life cycle consumption of coal by both types of biomass is 
smaller than the coal consumption of the agricultural 
phase. 

3.2.2 Least squares equivalent partitioned matrices 

Reducing the number of products of the original data to the 
same dimension by instituting the recovery and substitution of 
inerts for sand, one can determine the set of allocation weights 
consistent with OLS allocation for a given vector of 


emissions, b, by reconciling Eqs. 9 and 10. That is, one solves 
the equation 

«)A= (A T ) + b (17) 

for the allocation weights which appear in A t and b„. In 
Table 3, the allocation weights for selected interventions are 
compared to the weights for the traditional allocations, as 
derived from information in (Marvuglia et al. 2010, Tables 3 
and 4). As crude oil, natural gas, and brick (after assuming 
substitution of inerts for sand) were derived from single func¬ 
tion processes, they are omitted from the table. For OLS 
allocation, two other interventions are presented along with 
coal to illustrate that implied weights may differ for each of the 
distinct interventions, as if a different allocation principle has 
been chosen in each case. Note that in all outcomes, the 
weights assigned to a pair of products sum to 1; however, 
some of the weights assigned to joint products by ordinary 
least squares are clearly indefensible in the sense that some 
inputs may become outputs and while other products may be 
assigned a share of inputs greater than 100 % in the corre¬ 
sponding technology and interventions matrices. 

Both DLS and TLS derive their inventories from a 
perturbed version of Eq. 1. For each type of intervention, 
DLS and TLS will yield a scalar intensity associated with 
brick r bnck . Fixing this value still leaves five degrees of 
freedom in a solution of Eq. 8. One could find infinitely many 
intensity vectors in which the component ;- brick is equal to DLS 
or TLS outcome, perhaps even with all positive weights. For 
example, one possible family of TLS-equivalent partitioning 
allocations for coal could be obtained by setting (w wh i te day. 
wWd. w 0 ii ve ca ke) = (0.5, 0.5, 0.8) and searching the level curve 
/- h rick=250 MJ coal of the function shown in Fig. 3 for corre¬ 
sponding (w e i ec tricity, Wdi ese i oil) pairs. 

The set of all solutions to the underdetermined system 
of equations Eq. 8 is infinite, and each allocation, least 
squares techniques now included, merely selects one 
particular solution and interprets the coordinate of inter¬ 
est in the solution vector. The utility of any of these 
allocation methods is questionable if no further thought 
is given to the side conditions associated with them. 
When acquisition of additional data is simply impossi¬ 
ble, any statement as to the quality or correctness of a 
particular inventory after allocation must be tied to just 
the belief that the side conditions implied by the 
partitioning criterion are at least plausible. It is impera¬ 
tive that a decision maker should at least try to under¬ 
stand what those implications are! 

3.3 Discrepancy vectors 

The reported discrepancy vectors for economic and physical 
allocations (both invertible systems) were given along with 
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Table 3 Description of stated or 
implied weights assigned to each 
pair of coproducts. Weights for 
economic and physical allocation 
were as stated in (Marvuglia et al. 
2010, Tables 3 and 4). In OLS 
allocation, all pairs of coproduct 
weights sum to 1, but some 
weights are negative or exceed 1 


Unallocated Individual Stated or implied allocation weights 

process products - 




Economic 

Physical 

OLS MJ 

geotherm. 

energy 

OLS kg 
COD 

OLS, 

MJ coal 

Electricity and heat 

MJ of electricity 

0.518 

0.8 

-0.375 

0.206 

-0.289 


MJ of heat 

0.482 

0.2 

1.375 

0.794 

1.289 

Clay 

kg of white clay 

0.524 

0.5 

1.600 

-0.921 

0.516 


kg of red clay 

0.476 

0.5 

-16.180 

312.459 

0.262 

Sand and gravel 

kg of gravel 

0.577 

0.5 

17.180 

-311.459 

0.738 


kg of olive cake 

0.423 

0.5 

-0.600 

1.921 

0.484 

Supply of biomass 

kg of sand 

0.713 

0.8 

-5.515 

0.494 

0.742 


kg of straw 

0.287 

0.2 

6.515 

0.506 

0.258 

Oil derivatives 

MJ of diesel oil 

0.058 

0.5 

-0.004 

0.017 

-0.005 


MJ of fuel oil 

0.942 

0.5 

1.004 

0.983 

1.005 


discrepancy vectors from each of the three regression-based 
methods in (Marvuglia et al. 2010, Table 12). Observing that 
residuals among the regression-based techniques were of 
much greater magnitude than either of the traditional alloca¬ 
tion methods, the following conclusion was reached in 
(Marvuglia et al. 2010, p. 1028): 

“On the basis of the observation of the discrepancy vectors 
alone, the traditional allocation based solutions may there¬ 
fore appear to be preferred to the solutions obtained using 
the regression techniques. However, a complete judgment 
should also take into account the observation of the differ¬ 
ences among the inventory vectors since the inventory 


The discrepancies which were tabulated provide vital in¬ 
formation and merit more thorough examination, however. 

Comparison of the residuals of the traditional partitioning 
allocations to those of any of the least squares techniques does 
not make any sense, since the former results due to machine 
precision in the numerical solution of an invertible system of 
equations while the latter is due to the lack of a unique solution 
and the choice of an approximate solution of the overdeter¬ 
mined system Eq. 1 . The tabulated residuals ((Marvuglia et al. 
2010, Table 12)) associated with OLS, DLS, and TLS corre¬ 
spond to the quantities As 0 ls -/brick, A A DLS s DLS , and 
AH xls s xls - Af TLS , respectively. Beyond any residuals associ¬ 
ated with the functional unit in OLS and TLS, the residuals of 


vector is the nal goal of die LC1 phase.” the technology matrix in 

Table 4 Percentage change (100 * ZVJ DLS !\A\) in cells of A implied by AA DLS 

DLS and TLS must be carefully 


Electricity 
and heat 

Clay 

(red and white) 

Sand and 
gravel 

Crude 

oil 

Oil 

derivatives 

Natural 

gas 

Supply of 
biomass 

Bricks and 
inerts 

Functional 

unit 

MJ of electricity 


95.6 

14.1 

- 

-99.2 

- 

+ 


- 

MJ of heat 


- 

- 

253.5 

3.7 

3.0 

- 


+ 

kg of white clay 

- 


+ 

- 

- 

- 

+ 


- 

kg of red clay 

+ 


- 

+ 

+ 

+ 

- 


+ 

kg of recycled inerts 

+ 

- 

- 

+ 

+ 

+ 

- 


+ 

kg of sand 

- 

+ 


- 

- 

- 

+ 


- 

kg of gravel 

+ 

- 


+ 

+ 

+ 

- 


+ 

kg of olive cake 

- 

+ 

+ 

- 

- 

- 



- 

kg of straw 

+ 

- 

- 

+ 

+ 

+ 



+ 

MJ of crude oil 

- 

+ 

+ 

-18.3 

-2.7 

- 

+ 

+ 

- 

MJ of diesel oil 

+ 

-2.2 


26.3 





+ 

MJ of fuel oil 

- 

+ 

+ 

- 

-12.1 

- 

+ 


- 

MJ of natural gas 


- 

- 

+ 

+ 


- 


+ 

Ton of bricks 











Changes to structural zeros are indicated by +/-. Numeric values indicate a percentage change greater than +/- 2% of the original cell value. Bold face 
denotes a flow which has changed direction after revision 
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examined. Tables 4 and 5 examine these adjustments to the 
data. While small in the sense that the perturbations minimize 
the Frobenius norm of the correction matrix [A4|A f], the 
adjusted entries potentially: 

• change known zeros to small amounts of inputs or outputs, 

• indicate large changes for some element, perhaps even chang¬ 
ing sign so that inputs become outputs or vice versa, and 

• correspond to a very different technology matrix than the 
one intended. 

In Tables 4 and 5, the corrections matrix \AA\ A/| for DLS 
and TLS are expressed as percentage changes are given for 
non-zero elements of the original technology matrix. The two 
tables differ primarily in the bottom row and the functional 
unit column since TLS permits changes to the stated function¬ 
al unit as well as the technology matrix. In each table, a blank 
cell indicates a change smaller than 2 % (in absolute value) of 
the cell’s original value. The signs +/- indicate changes of 
structural zeros into small amounts of outputs or inputs, re¬ 
spectively. After perturbation, the brick production unit ap¬ 
pears to generate crude oil, changes in six cells were greater 
than 10 % of the magnitude of the original cell entry, and one 
cell entailed a 250 % change and a reversal of sign, indicating 
that, in the perturbed technology matrix A +A A , heat became 
an output of the crude oil process rather than an input used to 
do work and pump oil. This is another manifestation of the 
shortfall in system-wide energy inputs mentioned in 
Section 3.1. Consequently, the scale factors derived by DLS 


and TLS methods correspond to technology matrices with 
little resemblance to the original data. 

Note that the perturbations in A A also depend on the 
stated functional unit. Consequently, a researcher could per¬ 
form DLS or TLS on a given data matrix for two different 
functional units and derive rather different sets of corrections 
to the technology matrix. As stated in Marvuglia et al. (2010), 
it is possible that the recorded data were measured with some 
error, but surely the extent of the measurement errors should 
be independent of which functional unit is stated for study. 

In summary, Marvuglia et al. appeared to favor the TLS 
method on the grounds that it (1) produced (a few) results that 
were similar to physical allocation and (2) on the basis of 
asymptotic properties of TLS parameter estimates (i.e., strong 
consistency of TLS estimators). Regarding the first point, as 
discussed in the previous section, one should not necessarily 
expect agreement among inventories derived under different 
constraints. Among the five tabulated inventories, one can find 
agreement in some components of the inventories and radical 
disagreement in others, leaving the issue of which value(s) 
seem most reasonable up for debate. The argument for the 
second point simply appears to be that TLS will yield better 
results on bigger data sets. This reasoning is neither justified 
theoretically, nor does it resolve the allocation problem in any 
way. To act on the TLS scale factors, one must believe that the 
changes to the technology matrix A+A.A IIS are plausible. 
These perturbations may be quite significant as they could 
entail changes to structural zeros and perhaps changes of inputs 
into outputs or vice versa. In any case, the asymptotic properties 


Table 5 Percentage change (100 * A4 TLS l\A \) in cells of A implied by By assumption, TLS pennits small discrepancies in the stated functional 

AT tls . The same sign conventions used in Table 4 apply. DLS and TLS unit; small changes are also introduced in the final row to compensate for 

differ primarily by perturbations in the final row of the technology matrix. error in the stated functional unit 



Electricity 
and heat 

Clay (red and 
white) 

Sand and 
gravel 

Crude 

oil 

Oil 

derivatives 

Natural 

gas 

Supply of 
biomass 

Bricks and 
inerts 

Functional 

unit 

MJ of electricity 


95.6 

14.1 

- 

-99.2 

- 

+ 


- 

MJ of heat 


- 

- 

253.5 

3.7 

3.0 

- 


+ 

kg of white clay 

- 


+ 

- 

- 

- 

+ 


- 

kg of red clay 

+ 


- 

+ 

+ 

+ 

- 


+ 

kg of recycled 

+ 

- 

- 

+ 

+ 

+ 

- 


+ 

inerts 










kg of sand 

- 

+ 


- 

- 

- 

+ 


- 

kg of gravel 

+ 

- 


+ 

+ 

+ 

- 


+ 

kg of olive cake 

- 

+ 

+ 

- 

- 

- 

- 


- 

kg of straw 

+ 

- 

- 

+ 

+ 

+ 



+ 

MJ of cmde oil 

- 

+ 

+ 

-18.3 

-2.7 

- 

+ 

+ 

- 

MJ of diesel oil 

+ 

-2.2 


26.3 





+ 

MJ of fuel oil 

- 

+ 

+ 

- 

-12.1 

- 

+ 


- 

MJ of natural gas 


- 

- 

+ 

+ 


- 


+ 

Ton of bricks 

- 

+ 

+ 

- 

- 

- 

+ 


- 


Changes to structural zeros are indicated by +/-. Numeric values indicate a percentage change greater than +/- 2% of the original cell value. Bold face 
denotes a flow which has changed direction after revision 
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of s-tls do not guarantee that the discrepancies are small in any 
practical sense. The issue of allocation is one of rank deficiency 
or non-identifiability of the inventory of interest. In fact, apply¬ 
ing TLS on a data set of the size of the entire unallocated 
ecoinvent vl.3 database with 2,632 rows (products) and 2,471 
columns (processes) (Marvuglia et al. 2010, p. 1021) would not 
overcome the rank deficiency, requiring some 161 additional 
constraints, in any meaningful way. 


4 Conclusions 

In actuality, the fundamental problem of allocation is one of 
too few linear equations and too many unknowns. The com¬ 
ponents of the intensity matrix to be interpreted in the inven¬ 
tory step may not be uniquely identified in the given data. In 
such cases, the final inventory value will be determined by the 
additional infonnation brought to bear on a solution of Eq. 8. 

Upon reexamination of the brick production data, issues of 
data quality and the representativeness of included processes 
may explain some of the unusual disparities in the originally 
provided inventories. Each allocation by least squares tech¬ 
niques or by more traditional means of partitioning will yield a 
result irrespective of data quality. Consequently, it is important 
to define what constitutes a legitimate technology matrix and 
to identify and correct any gross errors before carrying out the 
remaining steps of the life cycle analysis. 

Of more general concern, an inventory produced by 
partitioning must necessarily bring some additional infonna¬ 
tion to bear on the problem. In that sense, the least squares 
techniques are no different than more traditional partitioning 
allocations. Inventories generated by least squares techniques 
are, in fact, equivalent to some partitioning of the technology 
and interventions matrices. We caution practitioners that the 
use of least squares techniques for the purposes of allocation is 
not a black-box solution to the problem of allocation and may 
not generate an output that is any more physically meaningful 
or actionable than already existing methods of partitioning. 

A square, invertible technology matrix obtained through 
allocation often marks the starting point for structural path 
analysis and many existing methods of uncertainty propaga¬ 
tion in LCA. Like the brick example, sensitivity of the inven¬ 
tory to allocation assumptions is often examined through 
application of several allocations to the same data. The math¬ 
ematical structure of the allocation problem pennits more 
general and thorough understanding. In particular, forthcom¬ 
ing work will provide ways to address the sensitivity of 
specific comparisons of products to allocation assumptions 
Hanes et al. (Forthcoming) and describe conditions for the 
existence of inverse-positive outcomes Cruze (2013) which 
may facilitate the exploration of the set of all allocations. 


To satisfactorily avoid partitioning, practitioners must ei¬ 
ther (1) acquire more data on similar processes involved in the 
life cycle of the target product or upstream inputs or (2) 
state additional assumptions about life cycle intensities of 
select products within the system more judiciously. As it 
stands, practitioners can follow prescribed recipes for alloca¬ 
tion without regard to the veracity of additional constraints, 
much to the detriment of the inventory and subsequent impact 
assessment and interpretation phases of any life cycle 
assessment. 
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